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ABSTRACT 

We present spectral energy distribution (SED) models for the edge-on classical T Tauri star HH30 IRS 
that indicate dust grains have grown to larger than 50/im within its circumstellar disk. The disk geometry 
and inclination are known from previous modeling of multiwavelength Hubble Space Telescope images 
and we use the SED (0.5/im < A < 3mm) to constrain the dust size distribution. Model spectra are 
shown for different circumstellar dust models: a standard ISM mixture and larger grain models. As 
compared to ISM grains, the larger dust grain models have a shallower wavelength dependent opacity: 
smaller at short wavelengths and larger at long wavelengths. Models with the larger dust grains provide a 
good match to the observed SED of HH30 IRS. Although the currently available SED is poorly sampled, 
we estimate L-i, « 0.2Lq, Mdisk ~ 1-5 x 10"'^ Mq, and a power-law with exponential cutoff dust grain 
size distribution. This model provides a good fit to the currently available data, but mid and far-IR 
observations are required to more tightly constrain the size distribution. The accretion luminosity in our 
models is Lace 0■2L^, corresponding to an accretion rate M ^ 4 x 10~^Moyr~^. Dust size distributions 
that are simple power-law extensions (i.e., no exponential cutoff) yield acceptable fits to the optical/near- 
IR but too much emission at mm wavelengths and require larger disk masses up to Mdisk 0.5Mq. Such 
a simple size distribution would not be expected in an environment such as the disk of HH30 IRS, (i.e., 
where coagulation and accretion processes are occurring in addition to grain shattering), particularly over 
such a large range in grain sizes. However, its ability to adequately characterize the grain populations 
may be determined from more complete observational sampling of the SED in the mid to far-IR. 

Subject headings: radiative transfer — scattering — accretion, accretion disks — ISM: dust, extinction 
— stars: pre-main-sequence — stars: individual (HH30 IRS) 



1. INTRODUCTION 

The HST WFPC2 images of HH30 IRS (Burrows et al. 
1996) revealed the characteristic scattered light pattern 
of an edge-on disk in a classical T Tauri system: two al- 
most parallel nebulae separated by an opaque dust lane 
(e.g., Whitney & Hartmann 1992). Detailed modeling 
of the WFPC2 images enabled Burrows et al. to con- 
strain the size, shape, and optical depth of the disk. In 
addition, their best model images were constructed using 
dust grains with a scattering phase function which is more 
forward throwing than typical interstellar medium (ISM) 
grains, indicative of grain growth within the disk (Hansen 
& Travis 1974). 

Wood et al. (1998) presented scattered light models of 
edge-on disks showing that the width of the dust lane de- 
creases towards longer wavelengths where the dust is less 
opaque (assuming a wavelength dependence for absorption 
and scattering of typical ISM dust grains) . They predicted 
that the HH30 IRS nebulosity should appear almost point- 
like at the K band. However, the new NICMOS observa- 
tions of Cotera et al. (2001) show a distinct dust lane even 
at the longest NICMOS wavelength. Cotera et al. used 
the NICMOS images to study the near-IR wavelength de- 
pendence of the circumstellar dust opacity. They modeled 
the NICMOS images using a dust size distribution that 



has a smaller near-IR opacity and a shallower wavelength 
dependence than ISM dust. Their modeling showed that 
the extinction through the disk to the central source is 
Av = 7900. The derivation of this extinction requires de- 
tailed radiation transfer modeling, since the HST images 
are dominated by scattered light so standard reddenning 
estimates are not applicable. 

The scattered light models of Burrows et al.. Wood 
et al. (1998), and Cotera et al. fit the width of the 
dust lane, from which they derive Mdisk ^ 10~'*Mq. As 
noted by Burrows et al., what is really determined from 
scattered light modeling is the equatorial optical depth 
in the disk midplane, Teq oc KAMdisk- The mass deter- 
mination is therefore very uncertain since it depends on 
the adopted grain size distribution and opacity. Scattered 
light models of the edge-on disk system HK Tauri/c derive 
Mdisk ~ 10"'' Mq assuming ISM dust (Stapelfeldt et al. 
1998). However, DAlessio, Calvet, & Hartmann (2001) 
have shown that the edge-on disk of HK Tau/c may be 
modeled with a massive disk, Mdisk = O.OGSM©, in which 
the maximum grain size is 1 meter. Their results also 
demonstrate the utility of long wavelength observations 
for distinguishing between a low mass disk with ISM grains 
and a massive disk with a population of much larger grains. 

It is difficult to draw firm conclusions on circumstellar 
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dust properties since the often unknown inclination and 
disk geometry lead to degeneracies in SED models. How- 
ever, for HH30 IRS scattered light modeling provides the 
disk geometry and inclination. Therefore the circumstellar 
grain properties may be studied through detailed model- 
ing of its SED which covers the range 0.55/im < A < 3mm 
(Stapelfeldt & Moneti 1999). This paper presents model 
SEDs for IIII30 IRS that highlight the differences between 
ISM and larger dust grain models. Our modeling adopts 
the disk structure derived from fitting HST images and 
calculates the radiative equilibrium SED from a star plus 
disk system. We use our new Monte Carlo radiative equi- 
librium code for the calculations and show SED models 
at different inclination angles ranging from almost edge- 
on (appropriate to HII30 IRS) to pole-on viewing of the 
star -|- disk system. 

Section 2 presents the HH30 IRS SED taken from the 
literature, §3 describes the radiation transfer, disk model, 
and dust opacity, and §4 presents SED models for a passive 
reprocessing disk and also the effects of including viscous 
accretion luminosity. §5 summarizes our results. 

2. THE SPECTRAL ENERGY DISTRIBUTION OF HH30 IRS 

Our models are constructed with the goal of reproduc- 
ing the currently available SED for HH30 IRS. The SED 
displayed in Fig. 1 is taken from Mundt & Fried (1983: V 
and /), Vrba, Rydgren, & Zak (1985: H, K, and an upper 
hmit at L), and Reipurth et al. (1993: 1.3mm). The SED 
measured by ISO between 4.5/xm and 90/xm and additional 
mm fluxes (Stapelfeldt & Padgett 2001) were presented by 
Stapelfeldt & Moneti (1999, Fig. 3). 

The ISO data were also presented by Brandner et al. 
(2000) and their reduction produced slightly larger fluxes 
than that of Stapelfeldt & Moneti (1999). Both reduc- 
tions are shown on this figure. The fluxes in Fig. 1 were 
not taken simultaneously and HH30 IRS is known to vary 
in the optical by up to 1.4mag at V, but this amplitude 
decreases towards longcir wavelengths (Wood et al. 2000). 
Consequently, this variability will not affect our dust size 
distribution modeling since the different dust models are 
discriminated at long wavelengths (A ^ 20/im) where any 
variability is expected to be small. 

Although the 60yum and 90/um data are upper hmits, 
the SED shows the characteristic double-peaked spectrum 
of an edge-on disk: the direct starlight is obscured by the 
disk, the near-IR peak arises from scattered light, and the 
far-IR peak is due to reprocessing of stellar radiation in 
the disk and possibly the liberation of disk accretion lumi- 
nosity (e.g., Sonnhalter, Preibisch, & Yorke 1995; Boss & 
Yorke 1996; Men'shchikov & Henning 1997; D'Alessio et 
al. 1999). 

The larger flux measurement at 1.3mm is from Reipurth 
et al. (1993) and was observed with a 23" beam. The 
mm data from Stapelfeldt & Moneti (1999) are from high 
spatial resolution mapping of the HH30 IRS disk using 
the Owens Valley Radio Observatory (OVRO). As such 
the OVRO data probes emission from the disk as opposed 
to the unresolved data of Reipurth et al. (1993). The 
larger flux measurement of Reipurth et al. indicates that 
there is material beyond the disk observed by HST and 
OVRO, but the large beam size does not spatially resolve 
this emission or determine whether all of it is really asso- 
ciated with IIH30 IRS (see also the 1.3mm map of Motte 



& Andre 2001, their Fig. Ih). In our modeling, we are 
specifically interested in the grain properties within the 
resolved circumstellar disk and therefore our models fit 
the mm data of Stapelfeldt & Moneti (1999). The HST 
images are fit with i?disk ^ 200AU (Cetera et al. 2001) 
which is comparable to the extent of the mm data we are 
modefing (Stapelfeldt & Padgett 2001) and we adopt this 
outer radius (see § 3.2). Both 1.3mm points are shown in 
Fig. 1, but in subsequent figures only the OVRO data are 
shown. 

The foreground extinction to IIII30 IRS is not well de- 
termined as the central source is not directly visible at 
optical and near-IR wavelengths. Cotera et al. estimated 
Ay <4: and Fig. 1 shows dereddened SEDs for 1 < Ay < 4 
using the dereddening formulae of Cardelli, Clayton, & 
Mathis (1988) with Ry = 3.1. Comparing the dereddened 
SED to models of edge-on disks suggests Ay < 2. 

3. MODELS 

We now outline the radiation transfer technique, disk 
geometry, luminosity sources, and circumstellar dust prop- 
erties that are required in our modeling. 

3.1. Radiative Equilibrium Calculation 

The radiation transfer calculations use the Monte Carlo 
radiative equilibrium and temperature correction tech- 
nique of Bjorkman & Wood (2001) which we have adapted 
to simulate T Tauri disks (Bjorkman, Wood, & Whit- 
ney 2001, in preparation). The ouput of our code is the 
disk temperature structure (due to heating by stellar pho- 
tons and accretion luminosity) and the emergent SED at 
a range of viewing angles. A calculation of the hydrostatic 
disk structure (e.g., D'Alessio et al. 1999) can be included 
in the Monte Carlo technique, but this would require an 
iterative scheme. At present we have not implemented 
such a scheme and instead perform the radiative equilib- 
rium calculation within a fixed disk geometry as we now 
describe. 

3.2. Disk Structure 

We calculate model SEDs by performing the radiation 
transfer and radiative equilibrium calculation within the 
disk structure derived from scattered light modeling of the 
HST images of HH30 IRS. The HST images have been 
modeled by Burrows et al. (1996), Wood et al. (1998), 
and Cotera et al. (2001) with a fiared disk density struc- 
ture, 

P = Po(^^yexp-^[z/h{z,)]\ (1) 

where is the radial coordinate in the disk midplane, 
po is the density at the stellar surface and the scale height 
increases with radius, h = ho (m/R^,)^ . For our SED mod- 
eling we adopt: (3 = 1.25, a — 2.25, ho = 0.017i?*, giv- 
ing /i[100AU] = 17AU. The (3 and a values are consistent 
with detailed hydrostatic structure calculations (D'Alessio 
et al. 1999). The inner radius, i?o = QR*, is set to be the 
dust destruction radius for T.^ = 3500K and dust subli- 
mation temperature = 1600K. The outer disk radius 
is fixed at iidisk = 200AU to match the spatial extent of 
the HST images and the mm observations of Stapelfeldt & 
Padgett (2001). Therefore our SED modeling and derived 
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disk mass and dust properties apply to the material within 
-Rdisk = 200 AU which forms the observed SED. More dust 
likely exists beyond 200AU (Reipurth et al. 1993; Motte 
& Andre 2001), but it is not detected in the high spa- 
tial resolution HST or mm images and therefore does not 
contribute to the observed SED we are modeling. 

To fit the NICMOS F160W image, Cotera et al. (2001) 
used a slightly different geometry (/3 and a) from what we 
described above and derived Mdisk = 6.7 x 1O~*M0. This 
scattered light modeling probes the mass of small parti- 
cles that dominate the opacity at short wavelengths. As 
such, modeling of optical and near-IR images is insensitive 
to large particles, rocks, and planetesimals, and therefore 
this mass is a lower limit. 

3.3. Energy Sources 

The energy input to the disk is from stellar photons 
and accretion luminosity liberated in the disk. Initially 
we consider a passive disk heated by a uniformly bright 
star. § 4.3 places limits on the luminosity arising from 
viscous accretion. As discussed in the previous section, 
the disk structure is fixed for our radiation transfer simu- 
lations. We then use a disk theory (Shakura & Sunyaev 
1973) to determine the accretion rate for a given disk mass 
and structure. The accretion rate and viscosity parameter, 
M and adisk, are related to the disk parameters by 



M = \/l8^ adisk Vc po hl/R^ , 



(2) 



where the critical velocity Vc = ^/GM^JR^. The flux 
due to viscous disk accretion, GM^M /2R^, is generated 
throughout the disk midplane region according to 



dE 3GM^M 
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(3) 



We adopt = 0.5Mq and vary adisk to place limits on 
the accretion luminosity for a given A/disk- 

The variability of HH30's scattered light nebula (Bur- 
rows et al. 1996; Stapelfeldt et al. 1998; Cotera et al. 
2001) and ground based VRI photometry (Wood et al. 
2000) indicates non-uniform ilhunination of the disk due 
to either stellar hotspots or shadowing due to large in- 
homogeneities in the inner disk. However, D'Alessio et al. 
(1998) showed that hot accretion rings do not significantly 
alter the SED, so we assume a uniformly bright star and 
do not include disk warping in our models. 

3.4. Circum,stella,r Dust Properties 

Burrows et al. fitted their HH30 IRS HST images using 
a dust scattering phase function that was more forward 
throwing than ISM dust grains and the near-IR modeling 
by Cotera et al. derived a shallower wavelength dependent 
opacity than that of ISM grains. Both these results are in- 
dicative of circumstellar dust grains that are larger (in a 
statistically significant sense, e.g., cross-section weighted 
mean size) than typical ISM grains. Grain growth within 
protoplanetary disks has been inferred from mm obser- 
vations and the slope of long wavelength SEDs — mm 
fluxes are larger and the SED slopes are shallower than 
predicted from ISM grains (e.g., Beckwith et al. 1990; 



Beckwith & Sargent 1991). In order to examine the pos- 
sibility of grain growth within the IIII30 disk as well as 
the potential characteristics thereof, we construct models 
with several different dust size distributions: a standard 
ISM mixture and several distributions with larger grains 
including both simple power-laws and power-laws with an 
exponential cutoff. 

For the ISM dust, we adopt the size (mass) distribution 
of Kim, Martin, & Hendry (1994; KMH) for the canonical 
diffuse interstellar sightline (i.e, i?y=3.1). In addition, we 
also reproduce the model of Cotera et al. (2001). The 
interested reader is referred to these two works for specific 
details of the models. 

For our "larger grain" dust models, we combine a 
"power-law with exponential cutoff" (e.g., KMH) for com- 
ponents of amorphous carbon and astronomical silicates 
with solar abundance constraints for carbon and silicon 
(C/H ~ 3.5 X lO-'' and Si/H ~ 3.6 x lO'^, respectively; 
Anders & Grevesse 1989; Crevesse & Noel 1993). The pri- 
mary motivation for introducniig an additional complex- 
ity onto a simple power-law (one with a single exponent, 
e.g. D'Alessio et al. 2001) is that it is not expected to 
apply over such a large range of grain sizes. This has 
been indicated in analyses of both ISM and circumstel- 
lar dust grains (e.g., KMH; Jones, Tielens, & HoUenbach 
1996; Witt, Smith, & Dwek 2001; Suttner & Yorke 2001). 
In addition, rather than joining two separate power-laws 
(i.e., Witt, Smith, & Dwek 2001), we adopt a distribution 
that provides a smooth transition between two different 
regimes and that can reproduce the general behavior of 
a large range of interstellar sight lines (e.g., KMH). The 
functional form of this size distribution is : 



n(o) da = Citt ^ exp (— [a/ac]') da . 



(4) 



The parameters p, Uc, q, which control the distribution 
shape, are adjusted to fit the wavelength dependence of the 
observations. Of course, a simple power-law corresponds 
to the case oi q = oo. Ci is set by requiring the grains to 
completely deplete a solar abundance of carbon and sili- 
con. The subscript "z" refers to the different components 
of the dust model: amorphous carbon and silicate (see Ta- 
ble 1). This function is essentially the "modified gamma 
distribution" found in atmospheric and planetary science 
literature. The end-points of the size integration, amin and 
cimax, are also parameters. Not surprisingly, models with 
finite values of q are insensitive to the exact value of Omax 
as long as it is >> Cc. 

The opacity of the dust grains is modelled under the 
assumption of two separate populations of homogeneous, 
spherical particles. In general, we use published dielec- 
tric fimctions for astronomical silicate (Weingartner & 
Draine 2001) and amorphous carbon (BE type, Zubko et 
al. 1996). However, our desired wavelength range extends 
longward of both sets of dielectric functions and slightly 
shortward of the amorphous carbon values (i.e., 100-500 
A not in published set). We use a subtractive Kramers- 
Kronig analyis (cf. Ahrenkiel 1971; routine kindly sup- 
plied by K. Snook), specifying the imaginary index of re- 
fraction in the regions of interest. This approach also re- 
quires specification of the real index of refraction at some 
reference point; we use a value of the real refractive in- 
dex at a wavelength adjacent to the missing range from 
the published data. For the long wavelength extension. 
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we use the approximation that imaginary index is propor- 
tional to (e.g., Pollack ct al. 1994): iV = -1 for 
the silicates and A'' = 1/6 for the amorphous carbon. The 
latter choice departs from the recommendation of Pollack 
ct al. for metals (A^ = 1), but is a much better match 
for the slope of the published data in the far-IR region. 
Fortunately, our model fits arc fairly insensitive to the ex- 
act choice of A''. For the short wavelength extension, we 
use the imaginary index of the BEl of Rouleau & Martin 
(1991) scaled to match smoothly with that of Zubko et al. 

The scattering properties are calculated using a com- 
bination of Mie series solution and geometrical optics al- 
gorithms. For the former, we employ the code of Wis- 
combe (1996) while the latter regime is handled with the 
code of Macke and collaborators {rt- ellipsoid, i.e., Macke & 
Mishchenko 1996). The transition point between the two 
algorithms is made at a size parameter of 20,000. Though 
this choice is somewhat arbitrary, at this point the forward 
scattering properties as calculated by each algorithm dif- 
fer at a level below 1% (the phase matrix elements are 
in agreement at the 5-10% level). This is also true at a 
size parameter of 10,000, but we chose the larger value 
in order to allow the phase functions in the visible and 
NIR to be calculated entirely by the Mie algorithms; this 
avoids introducing the small-scale (^ 10% outside the cen- 
tral diffraction peak) high-frequency oscillations produced 
by the geometrical optics algorithm. 

Figures 2 and 3 show the spectral shape of the opac- 
ity, albedo, and cosine asymmetry parameter (i.e., g) for 
the various dust models we consider. The size; distribution 
parameters are listed in Table 1. Although we use Mie 
theory for calculating the scattering properties, in the ra- 
diative equilibrium calculation we adopt a simple Heyney- 
Greenstein scattering phase function characterized by the 
single parameter, g (Heyney & Greenstein 1941). 

4. MODEL SEDS 

A tentative spectral type of MO for the illuminating 
source of HH30 IRS was assigned by Kenyon et al. (1998). 
However, because the central star is totally obscured by 
the disk at optical and ncar-IR wavelengths, the spectral 
type is very uncertain. Scattered light models indicate 
PS O.2L0, so we adopt = 3500K (Kurucz model at- 
mosphere, Kurucz 1994) and i?* = 1.2ir!0. The distance 
to HH30 IRS is taken to be 140pc. The following sections 
show the resulting SEDs for different dust models, start- 
ing with a purely passive reprocessing disk. For each dust 
model we adopt the disk size and structure described in 
§3.2 and determine the disk mass from fitting the NICMOS 
F160W image. We then calculate the model SED for com- 
parison with the data presented in figure 1. Figure 4 shows 
the NICMOS F160W image and scattered light model fits 
for the five dust types described in §3.4 and summarized 
in Table 2. The Ay values in this table are calculated dif- 
ferently than those presented by Cotera et al. (2001). We 
have simply integrated through the disk along a viewing 
angle oi i = 84°, whereas Cotera et al. measured how 
much light escapes in the Monte Carlo simulation to get 
an effective Ay- 

The width of the dust lane in the scattered light im- 
ages is sensitive not only to the disk mass (and therefore 
Teqand Ah), but also to the scattering phase function. If 
g were the same for all models then Ah would be the 



same also. However, as the phase function becomes more 
forward throwing for the larger dust grain models (i.e., g 
increases), a larger optical depth is required to obtain the 
same width of the dust lane. 

4.1. Reprocessing Disk SEDs for i = 84° 

Figures 5 and 6 show the HH30 IRS data and model 

SEDs (for i = 84° derived from scattered light models) for 
the dust types described in §3.4 and summarized in Ta- 
ble 1. The models show the characteristic double peaked 
SEDs of highly inclined disks. The disk mass has been 
chosen so that each model matches the NICMOS F160W 
image (Cotera ct al. 2001). The different dust models 
are distinguished from one another at longer wavelengths. 
Although the ISM grains and the Cotera et al. grains pro- 
duce similar mm fluxes, they differ in the 1 — 10/im range 
where the ISM grains have a steeper wavelength depen- 
dent opacity that yields too narrow a dust lane in NICMOS 
scattered light images. The mm fluxes are best reproduced 
with larger grain size distributions. Our pure reprocessing 
disk models yield a rather low flux in the 1 — 10//m range 
and §4.3 shows models including accretion luminosity that 
provide a better fit to this spectral regime. 

The different models, disk mass estimates, midplane op- 
tical depths (Teg), and extinction to the star for i = 84° 
are summarized in Tables 1 and 2. The simple power-law 
size distributions can have very large mm optical depths. 
While dust size distributions with small Omax can be used 
to fit individual scattered light images, the mm fluxes 
are best fit with the large grain models. These models 
require very different masses and our SED models show 
that they may be distinguished in the mid to far-IR mea- 
sured by SIRTF and in the sub-mm measured by SMA 
and ALMA. That is to say, although a simple power-law 
is not physically plausible in objects like HH30 IRS, the 
issue of whether it adequately reproduces the SED can be 
determined with sufficient observation coverage. 

4.2. Reprocessing Disk SEDs for All Inclinations 

Scattered light modeling of HH30 IRS indicates that 
we are seeing a Classical T Tauri star viewed almost edge- 
on. Whitney & Hartmann (1992) presented scattered light 
disk models for a range of viewing angles showing that 
edge-on disks, where the central star is obscured, will be 
most easily detectable. Our edge-on SED models agree 
well with the current HH30 IRS data and Fig. 7 shows 
model SEDs at a range of viewing angles. As expected, 
for viewing angles where the central source is miobscured, 
our HH30 IRS models do indeed resemble the SED of a 
Classical T Tauri star. The silicate features at 9/tim and 
18/im are most prominent in the ISM dust model. As 
the grain size increases and the opacity is dominated by 
larger grains, the relative strengths of the silicate features 
decrease (Fig. 2). This is reflected in the SEDs with the 
ttc = 50/xm size distribution showing a very smooth SED 
across the silicate features. However, at present our mod- 
els do not allow for different dust compositions at different 
locations within the disk as required for a detailed investi- 
gation of the disk chemistry and mineralogy through anal- 
ysis of crystalline features (e.g., Chiang et al. 2001). 

4.3. Effects of Accretion Luminosity 
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HH30 IRS is tlie source of a powerful jet (Mundt & 
Fried 1983; Mundt ct al. 1990; Graham & Hcycr 1990; 
Lopez et al. 1995), presumably driven by accretion lumi- 
nosity. The low disk mass for HH30 IRS derived using 
ISM dust models is unlikely to support a large accretion 
rate in the context of a disk theory. Therefore the larger 
disk masses we derive arc more appealing as they allow for 
significant accretion rates. Figure 8 shows models which 
include disk heating from viscous accretion for Model 1, 
which has Afdisk = 1-5 x IO'^Mq. 

We find iacc = 0.2L^ provides a good match to the 
1 — 10/um fluxes presented by Brandner et al. (2000). As- 
suming Mi, = 0.5Mq then M = 4 x lO'^^Moyr'^ and 
Q^disk = 0.02. A larger accretion luminosity produces too 
much 1 — 10/im emission compared to the observations. 
Figure 8 also shows models in which the accretion luminos- 
ity is lower (Lacc/J^* = 5%, 2%). Pole-on viewing yields 
typical near-IR excesses arising from accretion disks. 

The large uncertainties in the 1 — 10/im fluxes in addi- 
tion to the unknown and adisk preclude an accurate 
determination of the accretion rate. However, the accre- 
tion rate we derive is comparable to the total outflow rate 
in the jets, Mjct =■ 3.5 x lO^^MQyr^^ estimated by Bac- 
ciotti, Eisloffel, & Ray (1999). The combination of SED 
models of more accurate SIRTF fluxes and spectral line 
accretion diagnostics (e.g., GuUbring et al. 1998) are re- 
quired for a more accurate determination of the accretion 
rate. 

4.4. Components of the SED 

A feature of Monte Carlo radiation transfer is that we 
may track the history of each individual photon packet 
within the simulation and determine the relative contri- 
butions to the total SED of direct, scattered, reprocessed, 
and accretion luminosity. Figure 9 shows model SEDs for 
dust Model 1 at inclinations i = 0°, i = 60°, i = 71°, and 
i = 84° along with their component parts. The scattered 
light comprises stellar photons that have been (multiply) 
scattered; thermal photons are stellar photons that have 
been absorbed, reprocessed, and possibly scattered by the 
disk; accretion photons are those that initiated as accre- 
tion luminosity, but may have subsequently been scattered 
or absorbed and re-radiated. 

For pole-on viewing, the SED is dominated by direct 
starlight at optical wavelengths and thermal (i.e., repro- 
cessed) disk emission at long wavelengths. Accretion lumi- 
nosity contributes significantly at near and mid-IR wave- 
lengths, while scattered starlight contributes a few percent 
of the SED in the optical and near-IR and much less at 
longer wavelengths. Scattering will always yield an excess 
above the photospheric level (see also Whitney & Hart- 
mann 1992). For edge-on viewing, the situation is differ- 
ent: the star is totally obscured by the disk and scattered 
starlight dominates the SED in the optical and near-IR. 

5. SUMMARY 

We have presented Monte Carlo radiation transfer cal- 
culations for the SED of HH30 IRS. Our model com- 
prises a central stellar source surrounded by a circum- 
stellar disk whose structure was derived from modeling 
HST images. The dust opacity in the disk has a shallower 
wavelength dependence than typical ISM dust and the av- 



erage grain size is significantly larger. The model that 
matches the observed SED has a dust size distribution ex- 
tending to larger grain sizes than ISM dust; specifically 
cimax = 1mm with an exponential cutoff scale-length of 
flc = 50/im. The opacity at long wavelengths is roughly 
K ^ A^^, in agreement with that derived from sub- mm 
continuum SEDs in many classical T Tauri stars (Beck- 
with et al. 1990). While size distributions with most of 
the mass in smaller grains can be constructed to fit optical 
and near-IR images, these grain models produce too little 
mm emission. The disk mass we derive using the model 
above is Afdisk = 1.5 x 10^'^ M©, considerably larger than 
the Mdisk = 2 X lO-^M© found by Wood et al. (1998) 
using ISM grains. 

The derived disk mass is a lower limit, since the SED 
is not very sensitive to large particles, such as rocks and 
planetesimals. Estimations of disk mass from mm obser- 
vations are uncertain due to the poorly known opacity at 
mm wavelengths. For HH30's disk, Reipurth et al. (1993) 
estimated A'/disk = 0.02A'/q if Ki.smm = 0.003cm^/g or 
Mdisk = 0.003AfQ if Ki.3mm = 0.02cmVg- The opacity for 
our grain size distribution is Ki.smm = 0.08cm^/g, again 
demonstrating the sensitivity of the disk mass estimate to 
the adopted mm opacity. Observations with SIRTF and 
at sub-mm wavelengths will give a more complete picture 
of the SED of HH30 IRS and other edge-on disks, allow- 
ing us to further constrain the dust size distribution and 
ascertain the extent of the grain growth within the disk. 
In particular, the SED in the 20 — 800/im range will dis- 
criminate between the dust model with exponential cut-off 
(corresponding to Mdisk = 1-5 x 10~^Mq) and the sin- 
gle power-law size distributions which require much more 
massive disks (Afdisk ^ 0.5AfQ). 

Grain growth in protoplanetary disks has been in- 
ferred from observations: e.g., shallower slopes of the 
long wavelength SEDs and larger mm fiuxes compared 
with those expected for ISM dust parameters (Beckwith 
ct al. 1990; Beckwith & Sargent 1991); the observed wave- 
length dependence of dust lane widths in near-IR images 
of HH30 IRS (Cotera et al. 2001). Our SED modeling of 
HH30 IRS provides further evidence for grain growth and 
suggests larger grain sizes than those derived from scat- 
tered light models only. Grain growth within protoplan- 
etary disks is expected on timescales of 10^ years (Beck- 
with, Henning, & Nakagawa 2000 and references therein). 
The simulations of Suttner & Yorke (2001) find that com- 
pact grains can grow to a = 50/im within ^10^ years. In 
their model, however, the largest grains are found primar- 
ily in the disk midplane regions and as such will not have 
a big influence on the infrared SED, which is dominated 
by the upper layers or disk "atmosphere." Although the 
current implementation of our code adopts the same grain 
size distribution throughout the disk, our SED modeling 
is consistent with grain growth to larger than 50/im, even 
in the upper layers of the disk atmosphere that dominate 
the infrared SED. 

The mass that we derive for HH30's disk is more ap- 
pealing than the low disk mass derived from ISM dust 
models which are unlikely to support large accretion rates. 
More massive disks support larger accretion rates and for 
IIII30 IRS we infer an accretion luminosity of iacc ~ 
0.2i*. Our derived accretion rate, M = 4x lO~^Af0yr~^, 
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is uncertain as it depends on the unknown M^, and adisk 
viscosity parameter and the accurate determination of M 
requires spectral line observations in addition to SED mod- 
eling. 
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Table 1 
Dust Model Parameters 



Dust Model 


Component 




P 


q 


ac(y^m) 


amax(/tini) 


Cotera 


amc 


5 X IQ-^^ 


3.5 


1.0 


0.55 


20 




sil 


6 X 10"^^ 


3.0 


1.0 


0.55 


20 


Model 1 . . . . 


amc 


1.32 X 10-1^ 


3.0 


0.6 


50 


1000 




sil 


1.05 X 10-1^ 


3.0 


0.6 


50 


1000 


Model 2.... 


amc 


1.32 X 10-1^ 


3.5 


DO 




1000 




sil 


1.05 X 10"" 


3.5 


00 




1000 


Model 3 . . . . 


amc 


1.32 X 10"" 


3.0 


00 




1000 




sil 


1.05 X 10-1^ 


3.0 


oo 




1000 



The size distribution for the KMH model is taken from their mass distribution 
for the fit to the canonical (i?y=3.1) average diffuse ISM sight line. 



Table 2 
Disk Masses and Extinctions 



Dust Model 




Ah 




9{H) 


Av 


^1.3mm 


KMH 


3.5 X lO-'^ 


800 


31000 


0.28 


4600 


0.05 


Cotera 


3.5 X 10-4 


1000 


39000 


0.49 


2800 


0.05 


Model 1 . . . . 


1.5 X 10-3 


1800 


65000 


0.58 


3300 


7.9 


Model 2 . . . . 


0.015 


3700 


130000 


0.60 


6500 


120 


Model 3 . . . . 


0.5 


13000 


440000 


0.75 


15000 


4200 



Optical depths are in the disk midplane, extinctions are along i = 84°. 
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Fig. 1. — Spectral energy distribution of HH30 IRS. Triangles arc the ISO reductions of Brandner et al. (1999). The ISO fluxes at 60;um 
and 90/im axe upper limits. Also shown are dereddened SEDs (stars) assuming 1 < Ay < 4, with Ay = 1 the lowest curve. 
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Fig. 2. Dust parameters for ISM grains (dotted lines), size distribution used by Cotera et al. (dashed lines), and the size distribution 
(Model 1) we use for fitting the HH30 SED (solid lines). 
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Fig. 3. — Dust parameters for the exponential cutoff dust model (solid line, Model 1) and simple power-law models, p = 3.5 (dashed line. 
Model 2),p = 3.0 (dotted line, Model 3). 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



Fig. 4. — NICMOS F160W image of HH30 IRS and scattered light models for the different dust size distribution models investigated in 
this paper. Axes are labeled in arcsecs, where larcsec = 140AU. 
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1 10 100 1000 

Fig. 6. — Data and model SEDs for the three different dust models presented in Fig. 3. Stars show dereddened fluxes assuming Ay = 1. 
Model SEDs are for the exponential cutoff (solid line, Model 1), and simpl^Piwer-laws, p = 3.5 (dashed line, Model 2), p = 3.0 (dotted line. 
Model 3). The input stellar spectrum is also shown. The models may be discriminated by the SED in the range 20/im < A < 1mm. 
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1 10 100 1000 

A(/zm) 



Fig. 7. — Model SEDs at a range of viewing angles for KMH grains (upper panel), Cotera et al. (middle panel) and the size distribution 
we use for fitting the HH30 IRS SED (lower panel, Model 1). Notice the decrease in the strength of the silicate features at 9 and 18;um for 
the larger grain models. 




Fig. 8. — Data and model SEDs including accretion luminosity. The model SEDs are for i = 84° (lower curves) and pole-on models (upper 
curves). For the different viewing angles going from upper to lower curves the accretion luminosity is Lacc/^'* = 0.2, 0.05, 0.05, and 0.0. 
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A(,un) 

Fig. 9. — Model showing the component parts of the SED for dust Model 1 at inclinations i = 0°, i = 60°, i = 71°, and i = 84° (lower 
panel). Each panel shows the input stellar spectrum (short dash), total SED (uppermost solid line), direct photons (dots), scattered starlight 
(lowest solid line), photons reprocessed by the disk (long dash), and accretion luminosity (dash dot dot). Direct starlight is not detectable for 
i = 84° and scattered starlight dominates the SED at optical and near-IR wavelengths where the star is totally obscured by the flared disk. 
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CONSTRAINING THE CIRCUMSTELLAR DUST SIZE DISTRIBUTION 
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Submitted to the Astrophysical Journal 

ABSTRACT 

We present spectral energy distribution (SED) models for the edge-on classical T Tauri star HH30 IRS 
that indicate dust grains have grown to larger than 50^m within its circumstellar disk. The disk geometry 
and inclination are known from previous modeling of multiwavclength Hubble Space Telescope images 
and we use the SED (0.5/im < A < 3mm) to constrain the dust size distribution. Model spectra are 
shown for different circumstellar dust models: a standard ISM mixture and larger grain models. As 
compared to ISM grains, the larger dust grain models have a shallower wavelength dependent opacity: 
smaller at short wavelengths and larger at long wavelengths. Models with the larger dust grains provide a 
good match to the observed SED of HH30 IRS. Although the currently available SED is poorly sampled, 
we estimate « 0.2Lq, Mdisk ~ 1-5 x IO^^Mq, and a power-law with exponential cutoff dust grain 
size distribution. This model provides a good fit to the currently available data, but mid and far-IR 
observations are required to more tightly constrain the size distribution. The accretion luminosity in our 
models is Lace <, 0-2L^ corresponding to an accretion rate M <^4:X lO^^MQyr^^. Dust size distributions 
that are simple power-law extensions (i.e., no exponential cutoff) yield acceptable fits to the optical/near- 
IR but too much emission at mm wavelengths and require larger disk masses up to M^isk ~ O-^Mq. Such 
a simple size distribution would not be expected in an environment such as the disk of HH30 IRS, (i.e., 
where coagulation and accretion processes are occurring in addition to grain shattering) , particularly over 
such a large range in grain sizes. However, its ability to adequately characterize the grain populations 
may be determined from more complete observational sampHng of the SED in the mid to far-IR. 
Subject headings: radiative transfer — scattering — accretion, accretion disks — ISM: dust, extinction 
— stars: pre-main-sequence — stars: individual (HH30 IRS) 



1. INTRODUCTION 

The HST WFPC2 images of HH30 IRS (Burrows et al. 
1996) revealed the characteristic scattered light pattern 
of an edge-on disk in a classical T Tauri system: two al- 
most parallel nebulae separated by an opaque dust lane 
(e.g., Whitney & Hartmann 1992). Detailed modeling 
of the WFPC2 images enabled Burrows et al. to con- 
strain the size, shape, and optical depth of the disk. In 
addition, their best model images were constructed using 
dust grains with a scattering phase function which is more 
forward throwing than typical interstellar medium (ISM) 
grains, indicative of grain growth within the disk (Hansen 
& Travis 1974). 

Wood et al. (1998) presented scattered light models of 
edge-on disks showing that the width of the dust lane de- 
creases towards longer wavelengths where the dust is less 
opaque (assuming a wavelength dependence for absorption 
and scattering of typical ISM dust grains). They predicted 
that the HH30 IRS nebulosity should appear almost point- 
like at the K band. However, the new NICMOS observa- 
tions of Cotera et al. (2001) show a distinct dust lane even 
at the longest NICMOS wavelength. Cotera et al. used 
the NICMOS images to study the near-IR wavelength de- 
pendence of the circumstellar dust opacity. They modeled 
the NICMOS images using a dust size distribution that 



has a smaller near-IR opacity and a shallower wavelength 
dependence than ISM dust. Their modeling showed that 
the extinction through the disk to the central source is 
Av = 7900. The derivation of this extinction requires de- 
tailed radiation transfer modeling, since the HST images 
are dominated by scattered light so standard reddenning 
estimates are not applicable. 

The scattered light models of Burrows et al.. Wood 
et al. (1998), and Cotera et al. fit the width of the 
dust lane, from which they derive Mdisk ^ 10^'^ Mq. As 
noted by Burrows et al., what is really determined from 
scattered light modeling is the equatorial optical depth 
in the disk midplane, Teq oc KxMdisk- The mass deter- 
mination is therefore very uncertain since it depends on 
the adopted grain size distribution and opacity. Scattered 
light models of the edge-on disk system HK Tauri/c derive 
Mdisk ~ 10^'' Mq assuming ISM dust (Stapelfeldt et al. 
1998). However, DAlessio, Calvct, & Hartmann (2001) 
have shown that the edge-on disk of HK Tau/c may be 
modeled with a massive disk, Mdisk = O.O65M0, in which 
the maximum grain size is 1 meter. Their results also 
demonstrate the utility of long wavelength observations 
for distinguishing between a low mass disk with ISM grains 
and a massive disk with a population of much larger grains. 

It is difficult to draw firm conclusions on circumstellar 
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dust properties since the often unknown inclination and 
disk geometry lead to degeneracies in SED models. How- 
ever, for HH30 IRS scattered light modehng provides the 
disk geometry and inclination. Therefore the circumstellar 
grain properties may be studied through detailed model- 
ing of its SED which covers the range 0.55jUm < A < 3mm 
(Stapelfeldt & Moneti 1999). This paper presents model 
SEDs for HH30 IRS that highlight the differences between 
ISM and larger dust grain models. Our modehng adopts 
the disk structure derived from fitting HST images and 
calculates the radiative equilibrium SED from a star plus 
disk system. We use our new Monte Carlo radiative equi- 
hbrium code for the calculations and show SED models 
at different inclination angles ranging from almost edge- 
on (appropriate to IIII30 IRS) to pole-on viewing of the 
star + disk system. 

Section 2 presents the HH30 IRS SED taken from the 
literature, §3 describes the radiation transfer, disk model, 
and dust opacity, and §4 presents SED models for a passive 
reprocessing disk and also the effects of including viscous 
accretion luminosity. §5 summarizes our results. 

2. THE SPECTRAL EXERGY DISTRIBUTION OF IIII30 IRS 

Our models are constructed with the goal of reproduc- 
ing the currently available SED for HH30 IRS. The SED 
displayed in Fig. 1 is taken from Mundt & Fried (1983: V 
and /), Vrba, Rydgren, & Zak (1985: iJ, K, and an upper 
hmit at L), and Reipurth et al. (1993: 1.3mm). The SED 
measured by ISO between 4.5/im and 90/im and additional 
mm fluxes (Stapelfeldt & Padgett 2001) were presented by 
Stapelfeldt & Moneti (1999, Fig. 3). 

The ISO data were also presented by Brandner et al. 
(2000) and their reduction produced slightly larger fluxes 
than that of Stapelfeldt & Moneti (1999)." Both reduc- 
tions are shown on this figure. The fiuxes in Fig. 1 were 
not taken simultaneously and HH30 IRS is known to vary 
in the optical by up to 1.4mag at V , but this amplitude 
decreases towards longer wavelengths (Wood et al. 2000). 
Consequently, this variability will not affect our dust size 
distribution modehng since the different dust models are 
discriminated at long wavelengths (A ^ 20/im) where any 
variability is expected to be small. 

Although the 60/^m and 90/im data are upper limits, 
the SED shows the characteristic double-peaked spectrum 
of an edge-on disk: the direct starlight is obscured by the 
disk, the near-IR peak arises from scattered hght, and the 
far-IR peak is due to reprocessing of stellar radiation in 
the disk and possibly the liberation of disk accretion lumi- 
nosity (e.g., Sonnhalter, Preibisch, & Yorke 1995; Boss & 
Yorke 1996; Men'shchikov & Henning 1997; D'Alessio et 
al. 1999). 

The larger fiux measurement at 1.3mm is from Reipurth 
et al. (1993) and was observed with a 23" beam. The 
mm data from Stapelfeldt & Moneti (1999) are from high 
spatial resolution mapping of the HH30 IRS disk using 
the Owens VaUey Radio Observatory {OVRO). As such 
the OVRO data probes emission from the disk as opposed 
to the unresolved data of Reipurth et al. (1993). The 
larger flux measurement of Reipurth et al. indicates that 
there is material beyond the disk observed by HST and 
OVRO, but the large beam size does not spatially resolve 
this emission or determine whether all of it is reafly asso- 
ciated with HH30 IRS (see also the 1.3mm map of Motte 



& Andre 2001, their Fig. Ih). In our modeling, we are 
specifically interested in the grain properties within the 
resolved circumsteflar disk and therefore our models fit 
the mm data of Stapelfeldt & Moneti (1999). The HST 
images are fit with -Raisk ~ 200AU (Cotera et al. 2001) 
which is comparable to the extent of the mm data we are 
modeling (Stapelfeldt & Padgett 2001) and we adopt this 
outer radius (see § 3.2). Both 1.3mm points are shown in 
Fig. 1, but in subsequent figures only the OVRO data are 
shown. 

The foreground extinction to HH30 IRS is not well de- 
termined as the central source is not directly visible at 
optical and near-IR wavelengths. Cotera et al. estimated 
Ay < 4: and Fig. 1 shows dereddened SEDs for 1 < Ay < 4 
using the dereddening formulae of Cardelli, Clayton, & 
Mathis (1988) with Ry = 3.1. Comparing the dereddened 
SED to models of edge-on disks suggests Ay < 2. 

3. MODELS 

We now outline the radiation transfer technique, disk 
geometry, luminosity sources, and circumstellar dust prop- 
erties that are required in our modeling. 

3.1. Radiative Equilibrium Calculation 

The radiation transfer calculations use the Monte Carlo 
radiative equilibrium and temperature correction tech- 
nique of Bjorkman & Wood (2001) which we have adapted 
to simulate T Tauri disks (Bjorkman, Wood, & Whit- 
ney 2001, in preparation). The ouput of our code is the 
disk temperature structure (due to heating by stellar pho- 
tons and accretion luminosity) and the emergent SED at 
a range of viewing angles. A calculation of the hydrostatic 
disk structure (e.g., D'Alessio et al. 1999) can be included 
in the Monte Carlo technique, but this would require an 
iterative scheme. At present we have not implemented 
such a scheme and instead perform the radiative equilib- 
rium calculation within a fixed disk geometry as we now 
describe. 

3.2. Disk Structure 

We calculate model SEDs by performing the radiation 
transfer and radiative equihbrium calculation within the 
disk structure derived from scattered light modeling of the 
HST images of HH30 IRS. The HST images have been 
modeled by Burrows et al. (1996), Wood et al. (1998), 
and Cotera et al. (2001) with a flared disk density struc- 
ture, 

P = Po{^^ exp-i[^//i(ro)]\ (1) 

where zu is the radial coordinate in the disk midplane, 
po is the density at the steflar surface and the scale height 
increases with radius, h = ho {uj I R^)^ . For our SED mod- 
eling we adopt: (5 = 1.25, a = 2.25, ho = 0.017i?^, giv- 
ing /i[100AU] = 17AU. The /3 and a values are consistent 
with detailed hydrostatic structure calculations (D'Alessio 
et al. 1999). The inner radius, Rq = 6R*, is set to be the 
dust destruction radius for T^- = 3500K and dust subli- 
mation temperature Td = 1600K. The outer disk radius 
is flxed at i?disk = 200AU to match the spatial extent of 
the HST images and the mm observations of Stapelfeldt & 
Padgett (2001). Therefore our SED modehng and derived 
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disk mass and dust properties apply to the material within 
-Rdisk = 200 AU which forms the observed SED. More dust 
Hkely exists beyond 200AU (Reipurth et a. 1993; Motte 
& Andre 2001), but it is not detected in the high spa- 
tial resolution HST or mm images and therefore does not 
contribute to the observed SED we are modeling. 

To fit the NICMOS F160W image, Cotera ct al. (2001) 
used a slightly different geometry {(3 and a) from what we 
described above and derived Mdisk = 6.7 x 10"'' Mq. This 
scattered light modeling probes the mass of small parti- 
cles that dominate the opacity at short wavelengths. As 
such, modeling of optical and near-IR images is insensitive 
to large particles, rocks, and planetesimals, and therefore 
this mass is a lower limit. 

3.3. Energy Sources 

The energy input to the disk is from stellar photons 
and accretion luminosity liberated in the disk. Initially 
we consider a passive disk heated by a uniformly bright 
star. § 4.3 places limits on the luminosity arising from 
viscous accretion. As discussed in the previous section, 
the disk structure is fixed for our radiation transfer simu- 
lations. We then use a disk theory (Shakura & Sunyaev 
1973) to determine the accretion rate for a given disk mass 
and structure. The accretion rate and viscosity parameter, 
M and Qdisk, are related to the disk parameters by 



M = Vl8^ Qdisk Vc po hl/R^ 



(2) 



where the critical velocity = \/GMZ[RI- The flux 
due to viscous disk accretion, GM^^M /2Ri,, is generated 
throughout the disk midplane region according to 



dE 
dAdt 



(3) 



We adopt = 0.5M(;5 and vary aaisk to place limits on 
the accretion luminosity for a given Mdisk- 

The variability of HH30's scattered light nebula (Bur- 
rows et al. 1996; Stapelfeldt et al. 1998; Cotera et al. 
2001) and ground based VRI photometry (Wood et al. 
2000) indicates non-uniform illumination of the disk due 
to either stellar hotspots or shadowing due to large in- 
homogeneities in the inner disk. However, D'Alessio et al. 
(1998) showed that hot accretion rings do not significantly 
alter the SED, so we assume a uniformly bright star and 
do not include disk warping in our models. 

3.4. Circumstellar Dust Properties 

Burrows et al. fitted their HH30 IRS HST images using 
a dust scattering phase function that was more forward 
throwing than ISM dust grains and the near-IR modeling 
by Cotera et al. derived a shallower wavelength dependent 
opacity than that of ISM grains. Both these results are in- 
dicative of circumstellar dust grains that are larger (in a 
statistically significant sense, e.g., cross-section weighted 
mean size) than typical ISM grains. Grain growth within 
protoplanetary disks has been inferred from mm obser- 
vations and the slope of long wavelength SEDs — mm 
fiuxes are larger and the SED slopes are shallower than 
predicted from ISM grains (e.g., Beckwith et al. 1990; 



Beckwitli & Sargent 1991). In order to examine the pos- 
sibility of grain growth within the HH30 disk as well as 
the potential characteristics thereof, we construct models 
with several different dust size distributions: a standard 
ISM mixture and several distributions with larger grains 
including both simple power-laws and power-laws with an 
exponential cutoff. 

For the ISM dust, we adopt the size (mass) distribution 
of Kim, Martin, & Hendry (1994; KMH) for the canonical 
diffuse interstellar sightline (i.e, i?y=3.1). In addition, we 
also reproduce the model of Cotera et al. (2001). The 
interested reader is referred to these two works for specific 
details of the models. 

For our "larger grain" dust models, we combine a 
"power-law with exponential cutoff" (e.g., KMH) for com- 
ponents of amorphous carbon and astronomical silicates 
with solar abundance constraints for carbon and silicon 
(C/H ~ 3.5 X 10-4 and Si/H ~ 3.6 x lO^'^, respectively; 
Anders & Grevesse 1989; Grevesse & Noel 1993). The pri- 
mary motivation for introducing an additional complex- 
ity onto a simple power-law (one with a single exponent, 
e.g. D'Alessio et al. 2001) is that it is not expected to 
apply over such a large range of grain sizes. This has 
been indicated in analyses of both ISM and circumstel- 
lar dust grains (e.g., KMH; Jones, Tielens, & Hollenbach 
1996; Witt, Smith, & Dwek 2001; Suttner & Yorke 2001). 
In addition, rather than joining two separate power-laws 
(i.e., Witt, Smith, & Dwek 2001), we adopt a distribution 
that provides a smooth transition between two different 
regimes and that can reproduce the general behavior of 
a large range of interstellar sight lines (e.g., KMH). The 
functional form of this size distribution is : 



n{a) da = Cia ^ exp (— [a/oc]') da . 



(4) 



The parameters p, ac, q, which control the distribution 
shape, are adjusted to fit the wavelength dependence of the 
observations. Of course, a simple power-law corresponds 
to the case of g = oo. C; is set by requiring the grains to 
completely deplete a solar abundance of carbon and sili- 
con. The subscript "«" refers to the different components 
of the dust model: amorphous carbon and silicate (see Ta- 
ble 1). This function is essentially the "modified gamma 
distribution" found in atmospheric and planetary science 
literature. The end-points of the size integration, a^in and 
flmax, are also parameters. Not surprisingly, models with 
finite values of q are insensitive to the exact value of Omax 
as long as it is >> Oc. 

The opacity of the dust grains is modelled under the 
assumption of two separate populations of homogeneous, 
spherical particles. In general, we use published dielec- 
tric functions for astronomical silicate (Weingartner & 
Draine 2001) and amorphous carbon (BE type, Zubko et 
al. 1996). However, our desired wavelength range extends 
longward of both sets of dielectric functions and slightly 
shortward of the amorphous carbon values (i.e., 100-500 
A not in pubUshed set). We use a subtractive Kramers- 
Kronig analyis (cf. Ahrenkiel 1971; routine kindly sup- 
plied by K. Snook), specifying the imaginary index of re- 
fraction in the regions of interest. This approach also re- 
quires specification of the real index of refraction at some 
reference point; we use a value of the real refractive in- 
dex at a wavelength adjacent to the missing range from 
the published data. For the long wavelength extension. 
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wc use the approximation that imaginary index is propor- 
tional to A"^' (e.g., Pollack et al. 1994): N = -1 for 
the sihcates and TV = 1/6 for the amorphous carbon. The 
latter choice departs from the recommendation of Pollack 
ct al. for metals {N = 1), but is a much better match 
for the slope of the published data in the far-IR region. 
Fortunately, our model fits are fairly insensitive to the ex- 
act choice of N. For the short wavelength extension, we 
use the imaginary index of the BEl of Rouleau & Martin 
(1991) scaled to match smoothly with that of Zubko ct al. 

The scattering properties are calculated using a com- 
bination of Mie series solution and geometrical optics al- 
gorithms. For the former, we employ the code of Wis- 
combe (1996) while the latter regime is handled with the 
code of Macke and collaborators {rt-ellipsoid, i.e., Macke & 
Mishchenko 1996). The transition point between the two 
algorithms is made at a size parameter of 20,000. Though 
this choice is somewhat arbitrary, at this point the forward 
scattering properties as calculated by each algorithm dif- 
fer at a level below 1% (the phase matrix elements are 
in agreement at the 5-10% level). This is also true at a 
size parameter of 10,000, but we chose the larger value 
in order to allow the phase functions in the visible and 
NIR to be calculated entirely by the Mie algorithms; this 
avoids introducing the small-scale (~ 10% outside the cen- 
tral diffraction peak) high-frequency oscillations produced 
by the geometrical optics algorithm. 

Figures 2 and 3 show the spectral shape of the opac- 
ity, albedo, and cosine asymmetry parameter (i.e., g) for 
the various dust models we consider. The size distribution 
parameters are hsted in Table 1. Although we use Mie 
theory for calculating the scattering properties, in the ra- 
diative equilibrium calculation we adopt a simple Heyney- 
Greenstein scattering phase function characterized by the 
single parameter, g (Heyney & Greenstein 1941). 

4. MODEL SEDS 

A tentative spectral type of MO for the ihuminating 
source of HH30 IRS was assigned by Kenyon et al. (1998). 
However, because the central star is totally obscured by 
the disk at optical and near-IR wavelengths, the spectral 
type is very uncertain. Scattered light models indicate 
« O.2L0, so we adopt = 3500K (Kurucz model at- 
mosphere, Kurucz 1994) and i?* = 1.2i?0. The distance 
to HH30 IRS is taken to be 140pc. The fohowing sections 
show the resulting SEDs for different dust models, start- 
ing with a purely passive reprocessing disk. For each dust 
model we adopt the disk size and structure described in 
§3.2 and determine the disk mass from fitting the NICMOS 
F160W image. We then calculate the model SED for com- 
parison with the data presented in figure 1. Figure 4 shows 
the NICMOS F160W image and scattered fight model fits 
for the five dust types described in §3.4 and summarized 
in Table 2. The Ay values in this table are calculated dif- 
ferently than those presented by Cotera et al. (2001). We 
have simply integrated through the disk along a viewing 
angle oi i = 84°, whereas Cotera et al. measured how 
much light escapes in the Monte Carlo simulation to get 
an effective Ay- 

The width of the dust lane in the scattered light im- 
ages is sensitive not only to the disk mass (and therefore 
Te^and Ah), but also to the scattering phase function. If 
g were the same for all models then Ah would be the 



same also. However, as the phase function becomes more 
forward throwing for the larger dust grain models (i.e., g 
increases), a larger optical depth is required to obtain the 
same width of the dust lane. 

4.1. Reprocessing Disk SEDs for i = 84° 

Figures 5 and 6 show the HH30 IRS data and model 
SEDs (for i = 84° derived from scattered light models) for 
the dust types described in §3.4 and summarized in Ta- 
ble 1. The models show the characteristic double peaked 
SEDs of highly inclined disks. The disk mass has been 
chosen so that each model matches the NICMOS F160W 
image (Cotera et al. 2001). The different dust models 
are distinguished from one another at longer wavelengths. 
Although the ISM grains and the Cotera et al. grains pro- 
duce similar mm fiuxes, they differ in the 1 — 10//m range 
where the ISM grains have a steeper wavelength depen- 
dent opacity that yields too narrow a dust lane in NICMOS 
scattered fight images. The mm fluxes are best reproduced 
with larger grain size distributions. Our pure reprocessing 
disk models yield a rather low flux in the 1 — 10/um range 
and §4.3 shows models including accretion luminosity that 
provide a better fit to this spectral regime. 

The different models, disk mass estimates, midplane op- 
tical depths (Teq), and extinction to the star for i = 84° 
are summarized in Tables 1 and 2. The simple power-law 
size distributions can have very large mm optical depths. 
While dust size distributions with small Cmax can be used 
to fit individual scattered light images, the mm fiuxes 
are best fit with the large grain models. These models 
require very different masses and our SED models show 
that they may be distinguished in the mid to far-IR mea- 
sured by SIRTF and in the sub-mm measured by SMA 
and ALMA. That is to say, although a simple power-law 
is not physically plausible in objects like HH30 IRS, the 
issue of whether it adequately reproduces the SED can be 
determined with sufficient observation coverage. 

4.2. Reprocessing Disk SEDs for All Inclinations 

Scattered light modeling of HH30 IRS indicates that 
we are seeing a Classical T Tauri star viewed almost edge- 
on. Whitney & Hartmann (1992) presented scattered light 
disk models for a range of viewing angles showing that 
edge-on disks, where the central star is obscured, wfil be 
most easily detectable. Our edge-on SED models agree 
well with the current HH30 IRS data and Fig. 7 shows 
model SEDs at a range of viewing angles. As expected, 
for viewing angles where the central source is unobscured, 
our HH30 IRS models do indeed resemble the SED of a 
Classical T Tauri star. The sificate features at 9/im and 
18/im are most prominent in the ISM dust model. As 
the grain size increases and the opacity is dominated by 
larger grains, the relative strengths of the silicate features 
decrease (Fig. 2). This is reflected in the SEDs with the 
flc = 50/im size distribution showing a very smooth SED 
across the silicate features. However, at present our mod- 
els do not allow for different dust compositions at different 
locations within the disk as required for a detailed investi- 
gation of the disk chemistry and mineralogy through anal- 
ysis of crystalfine features (e.g., Chiang et al. 2001). 

4.3. Effects of Accretion Luminosity 
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HH30 IRS is the source of a powerful jet (Mundt & 
Fried 1983; Mundt et al. 1990; Graham & Heyer 1990; 
Lopez et al. 1995), presumably driven by accretion lumi- 
nosity. The low disk mass for HH30 IRS derived using 
ISM dust models is unlikely to support a large accretion 
rate in the context of a disk theory. Therefore the larger 
disk masses we derive are more appealing as they allow for 
significant accretion rates. Figure 8 shows models which 
include disk heating from viscous accretion for Model 1, 
which has Afjisk = 1.5 x IO^^A/q. 

We find Z/acc = 0.2L^ provides a good match to the 
1 — 10/im fluxes presented by Brandner et al. (2000). As- 
suming = O.5M0 then M = 4 x lO'^^Afoyr"! and 
Q!disk = 0.02. A larger accretion luminosity produces too 
much 1 — 10/im emission compared to the observations. 
Figure 8 also shows models in which the accretion luminos- 
ity is lower (Lacc/-C'* = 5%, 2%). Pole-on viewing yields 
typical near-IR excesses arising from accretion disks. 

The large uncertainties in the 1 — 10/um fluxes in addi- 
tion to the unknown M^. and adisk preclude an accurate 
determination of the accretion rate. However, the accre- 
tion rate we derive is comparable to the total outflow rate 
in the jets, Mjet = 3.5 x lO~^M0yr~^ estimated by Bac- 
ciotti, Eisloffel, & Ray (1999). The combination of SED 
models of more accurate SIRTF fluxes and spectral line 
accretion diagnostics (e.g., Gullbring et al. 1998) are re- 
quired for a more accurate determination of the accretion 
rate. 

4.4. Components of the SED 

A feature of Monte Carlo radiation transfer is that we 
may track the history of each individual photon packet 
within the simulation and determine the relative contri- 
butions to the total SED of direct, scattered, reprocessed, 
and accretion luminosity. Figure 9 shows model SEDs for 
dust Model 1 at inclinations i = 0°, i = 60°, i = 71°, and 
i = 84° along with their component parts. The scattered 
light comprises stellar photons that have been (multiply) 
scattered; thermal photons are stellar photons that have 
been absorbed, reprocessed, and possibly scattered by the 
disk; accretion photons are those that initiated as accre- 
tion luminosity, but may have subsequently been scattered 
or absorbed and re-radiated. 

For pole-on viewing, the SED is dominated by direct 
starhght at optical wavelengths and thermal (i.e., repro- 
cessed) disk emission at long wavelengths. Accretion lumi- 
nosity contributes signiflcantly at near and mid-IR wave- 
lengths, while scattered starlight contributes a few percent 
of the SED in the optical and near-IR and much less at 
longer wavelengths. Scattering will always yield an excess 
above the photospheric level (see also Whitney & Hart- 
mann 1992). For edge-on viewing, the situation is differ- 
ent: the star is totally obscured by the disk and scattered 
starlight dominates the SED in the optical and near-IR. 

5. SUMMARY 

We have presented Monte Carlo radiation transfer cal- 
culations for the SED of HH30 IRS. Our model com- 
prises a central stellar source surrounded by a circum- 
stellar disk whose structure was derived from modeling 
HST images. The dust opacity in the disk has a shallower 
wavelength dependence than typical ISM dust and the av- 



erage grain size is signiflcantly larger. The model that 
matches the observed SED has a dust size distribution ex- 
tending to larger grain sizes than ISM dust; specifically 
flmax = 1mm with an exponential cutoff scale-length of 
Oc = 50jUm. The opacity at long wavelengths is roughly 
K ~ A^^, in agreement with that derived from sub-mm 
continuum SEDs in many classical T Tauri stars (Beck- 
with et al. 1990). While size distributions with most of 
the mass in smaller grains can be constructed to fit optical 
and near-IR images, these grain models produce too little 
mm emission. The disk mass we derive using the model 
above is Mjisk = 1.5 x IO^^Mq, considerably larger than 
the Mdisk = 2 X 10-*Mq found by Wood et al. (1998) 
using ISM grains. 

The derived disk mass is a lower limit, since the SED 
is not very sensitive to large particles, such as rocks and 
planetesimals. Estimations of disk mass from mm obser- 
vations are uncertain due to the poorly known opacity at 
mm wavelengths. For HH30's disk, Reipurth et al. (1993) 
estimated Mdisk = O.O2M0 if Ki..3nim = 0.003cm-^/g or 
Mdisk = O.OO3M0 if Ki.3mm = 0.02cmVg- The opacity for 
our grain size distribution is Ki.smm = 0.08cm^/g, again 
demonstrating the sensitivity of the disk mass estimate to 
the adopted mm opacity. Observations with SIRTF and 
at sub-mm wavelengths will give a more complete picture 
of the SED of HH30 IRS and other edge-on disks, allow- 
ing us to further constrain the dust size distribution and 
ascertain the extent of the grain growth within the disk. 
In particular, the SED in the 20 — 800//m range will dis- 
criminate between the dust model with exponential cut-off 
(corresponding to Mdisk = 1.5 x 1O~^M0) and the sin- 
gle power-law size distributions which require much more 
massive disks (Mdisk ^ 0.5A4p;). 

Grain growth in protoplanetary disks has been in- 
ferred from observations: e.g., shallower slopes of the 
long wavelength SEDs and larger mm fiuxes compared 
with those expected for ISM dust parameters (Beckwith 
et al. 1990; Beckwith & Sargent 1991); the observed wave- 
length dependence of dust lane widths in near-IR images 
of HH30 IRS (Cotera et al. 2001). Our SED modehng of 
HH30 IRS provides further evidence for grain growth and 
suggests larger grain sizes than those derived from scat- 
tered light models only. Grain growth within protoplan- 
etary disks is expected on timescales of 10^ years (Beck- 
with, Henning, & Nakagawa 2000 and references therein). 
The simulations of Suttner & Yorke (2001) find that com- 
pact grains can grow to a = 50//m within ~ 10^ years. In 
their model, however, the largest grains are found primar- 
ily in the disk midplane regions and as such will not have 
a big infiuence on the infrared SED, which is dominated 
by the upper layers or disk "atmosphere." Although the 
current implementation of our code adopts the same grain 
size distribution throughout the disk, our SED modeling 
is consistent with grain growth to larger than 50/um, even 
in the upper layers of the disk atmosphere that dominate 
the infrared SED. 

The mass that we derive for HH30's disk is more ap- 
pealing than the low disk mass derived from ISM dust 
models which are unlikely to support large accretion rates. 
More massive disks support larger accretion rates and for 
HH30 IRS we infer an accretion luminosity of L^cc ~ 
0.2L*. Our derived accretion rate, M = 4 x lO^^M0yr^^, 
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is uncertain as it depends on the unknown Af^ and ajisk 
viscosity parameter and the accurate determination of M 
requires spectral line observations in addition to SED mod- 
eling. 
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Table 1 
Dust Model Parameters 



Dust Model 


Component 




P 


q 




amax(/im) 


Cotera 


amc 


5 X 10-^^ 


3.5 


1.0 


0.55 


20 




sil 


6 X 10-18 


3.0 


1.0 


0.55 


20 


Model 1 . . . . 


amc 


1.32 X 10-1'^ 


3.0 


0.6 


50 


1000 




sil 


1.05 X 10-" 


3.0 


0.6 


50 


1000 


Model 2 . . . . 


amc 


1.32 X 10-1'^ 


3.5 


00 




1000 




sil 


1.05 X 10-1^ 


3.5 


00 




1000 


Model 3 . . . . 


amc 


1.32 X 10-1'^ 


3.0 


00 




1000 




sil 


1.05 X 10-1^ 


3.0 


00 




1000 



The size distribution for the KMH model is taken from their mass distribution 
for the fit to the canonical (J?y=3.1) average diffuse ISM sight Hne. 



Table 2 
Disk Masses and Extinctions 



Dust Model 




Ah 




9{H) 


Av 


^l.Smm 


KMH 


3.5 X 10"" 


800 


31000 


0.28 


4600 


0.05 


Cotera 


3.5 X 10-4 


1000 


39000 


0.49 


2800 


0.05 


Model 1 . . . . 


1.5 X 10-3 


1800 


65000 


0.58 


3300 


7.9 


Model 2 . . . . 


0.015 


3700 


130000 


0.60 


6500 


120 


Model 3 . . . . 


0.5 


13000 


440000 


0.75 


15000 


4200 



Optical depths are in the disk midplane, extinctions are along i = 84°. 
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Fig. 1. — Spectral energy distribution of HH30 IRS. Triangles are the ISO reductions of Brandner et al. (1999). The ISO fluxes at 60/im 
and 90/jm are upper limits. Also shown are dereddened SEDs (stars) assuming 1 < Ay < 4, with Ay ~ 1 the lowest curve. 
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Fig. 2. — Dust parameters for ISM grains (dotted lines), size distribution used by Cotera et al. (dashed lines), and the size distribution 
(Model 1) we use for fitting the HH30 SED (solid lines). 
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Fig. 3. — Dust parameters for the exponential cutoff dust model (solid line, Model 1) and simple power-law models, p = 3.5 (dashed line. 
Model 2), p = 3.0 (dotted line. Model 3). 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



Fig. 4. — NICMOS F160W image of HH30 IRS and scattered light models for the diflFerent dust size distribution models investigated in 
this paper. Axes are labeled in arcsecs, where larcsec — 140AU. 
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Fig. 7. — Model SEDs at a range of viewing angles for KMH grains (upper panel), Cotera et al. (middle panel) and the size distribution 
we use for fitting the HH30 IRS SED (lower panel, Model 1). Notice the decrease in the strength of the silicate features at 9 and for 
the larger grain models. 
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Fig. 8. — Data and model SEDs including accretion luminosity. The model SEDs are for i = 84° (lower curves) and pole-on models (upper 
curves). For the different viewing angles going from upper to lower curves the accretion luminosity is Lacc/^/* = 0.2, 0.05, 0.05, and 0.0. 
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Fig. 9. — Model showing the component parts of the SED for dust Model 1 at inclinations i = 0° , i — 60°, i = 71°, and i = 84° (lower 
panel). Each panel shows the input stellar spectrum (short dash), total SED (uppermost solid line), direct photons (dots), scattered starlight 
(lowest solid line), photons reprocessed by the disk (long dash), and accretion luminosity (dash dot dot). Direct starlight is not detectable for 
i — 84° and scattered starlight dominates the SED at optical and near-IR wavelengths where the star is totally obscured by the flared disk. 



